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Ames Research Center 


SUMMARY 


The possible occurrence of a viscous region near a surface from which 
there is rapid efflux of gas, accompanied by large heat transfer, is postu- 
lated and investigated theoretically. Such a viscous region, denoted as a 
boundary shock wave, may occur in the case of a large high-speed meteor enter- 
ing the earth T s atmosphere. The conditions across a boundary shock wave and 
its structure are calculated from the appropriate macroscopic equations 
reduced to closed-form expressions under the restriction of a perfect gas 
flowing at constant total enthalpy. 


INTRODUCTION 


The purposes of this report are to introduce the concept and to postulate 
the occurrence of boundary shock waves. A boundary shock wave may be defined 
as a thin region of viscous flow adjacent to a surface from which gas is flow- 
ing at a very high rate with large heat transfer. The flow in this region is 
dominated by heat conduction and viscous compressive stresses whose rapid 
decay causes rapid variations of the flow properties over a small distance. 

A detailed discussion of the causes of boundary shock waves is given in the 
section following the Introduction. 

If a boundary shock wave can occur as postulated, the distinct region of 
viscous flow would be possible only with high influx Reynolds number and with 
a source of high heat transfer (such as intense radiation being absorbed) at 
the wall. This high Reynolds number theory is the opposite limiting case to 
the more usual mass injection or vapor ablation at a rate small enough that 
the influx does not appreciably alter the boundary- layer character of an 
external flow along the wall. In the latter case of mass influx at a small 
rate (and small influx Reynolds number), the injected flow is entirely viscous, 
whereas injected flow at high influx Reynolds number must be essentially invis- 
cid a short distance from the surface. The boundary shock wave is therefore a 
"quick-transition region" belonging to the class of asymptotic phenomena dis- 
cussed by Friedrichs (ref. 1) . The boundary shock is analogous both to a 
boundary layer, in which the flow is along a wall, and to a shock wave, in 
which the flow is normal to the thin viscous region (see fig. 1; notation is 
defined in appendix A) . 

The boundary shock phenomenon may be important in understanding the flow 
over a certain class of meteoric fireballs . Large stone meteors that enter 
the earth 1 s atmosphere at very high velocity have a high rate of mass transfer 
from the surface by vapor ablation as a result of the intense radiative heat- 
ing of the surface by the extremely high temperature air behind the shock wave 
(see refs. 2 to 5 ) • When the ablation rate is so large that the distance, L, 



from the forward molten surface of the meteor to the interface between the air 
and the vaporized meteoric material is much greater than the "viscous length, " 

then Re^ = L/z v is large and a boundary shock wave is postulated to be 
present in the vapor flow, with inviscid flow just outside it. (See fig. 2, 
which represents the flow near the stagnation point of a blunt body.) 

A boundary shock wave may also occur in the flow of a gas out of a porous 
wall if there is a strong source of heat, such as radiation being absorbed at 
the surface and conducted back into the wall and gas in the pores . The gen- 
eral theory to be developed below will apply to that case if the flow can be 
laminar, not turbulent . 

The boundary shock wave will be analyzed here from the macroscopic, or 
continuum, point of view. The thermodynamic and flow variables must satisfy 
the classical macroscopic conservation laws, which will be approximated by 
using the Navier-Stokes equations of one- dimensional flow. The results in the 
general theory will be limited by the applicability of the approximate linear, 
i.e., first order, transport laws (the stress-strain rate relationship and 
the heat-conduction law) of the Navier-Stokes formulation. If the gradients 
become too steep, so that the gas is too far out of mechanical equilibrium, 
these transport laws will not give a very good approximation, and this must be 
kept in mind. It would be desirable to have a comparable microscopic analysis 
for comparison, but the state of the art in nonequilibrium kinetic theory for 
flow near surfaces may not at present be sophisticated enough to permit such 
an analysis . The one-dimensional macroscopic equations to be used here to 
study the flow through the boundary shock wave are the same as those that have 
been used successfully In the study of shock-wave structure. Professor 
Liepmann and his coworkers (ref. 6) have shown that the structure of most 
shock waves, especially in the downstream portion, is described very well by 
the Navier-Stokes equations. The structure of a boundary shock wave corre- 
sponds in many cases to a downstream portion of a simple-shock- wave solution. 
(Note then that the boundary shock may be entirely subsonic.) 

Calculation of the flow through a boundary shock will be based largely on 
the wealth of existing literature (see refs, 6 to 27) on shock waves and their 
structure. The relationships between the states and flow variables on both 
sides of a shock wave have analogous counterparts in relationships across a 
boundary shock which depend on the heat-transfer coefficient at the boundary 
as well as on the thermodynamic states (to be shown) . 

The sections which follow in this report will include: (l) a discussion 

of the causes of boundary shock waves; (2) development of the general theory 
for one-dimensional, laminar, steady flow through a boundary shock in a per- 
fect gas with longitudinal Prandtl number equal to unity, including explicit 
expressions for the jump conditions, the heat-conduction coefficient at the 
wall, and the structure, or variation of properties through the boundary shock; 
and ( 3 ) an approximate special application of the general theory to the spe- 
cial case of very rapid vapor ablation from a blunt body, with the properties 
of a stone meteor, entering the earth's atmosphere at high speed. 
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CAUSES OF BOimrARY SHOCK WAVES 


Consider a gas flowing out of a nearly .plane wall (fig. 1(c)) by either 
blowing through pores in the wall or vapor ablation, with radiative heat flux, 
-q r , being absorbed at the wall (to the left of b) in both cases. Let the 
influx Reynolds number, R%, be large enough that most of the flow within 
length L is inviscid. Since there is no source of viscous effects to the 
right of b, any viscous effects must originate at the wall (to the left of b). 
To understand the origin of the viscous effects in b to e, we must examine 
the mechanism of mass transfer within the region from s to b . 

Consider first the case of blowing through very small pores. At b the 
velocity has become essentially perpendicular to the wall. Assume that -q r 
is absorbed within a very short distance inside the surface s and is con- 
ducted back through the wall and into and back through the gas in the holes 
(to the left in fig. 1(c)) at a high rate. In a one- dimensional flow the 
viscous-compressive stress and heat- conduct ion flux are, respectively, 

t = p du/dx and q c = -k dT/dx (l) 

so that 

r/q, c = -Pr(&u/c p dT) (2) 

where 

Vv = qCp/k (3) 

and 

Vl = A + 2q = (V3V + (1/3 )k (4) 

and where Pr is the "longitudinal Prandtl number," q is the shear- viscosity 
coefficient, A is the second viscosity coefficient, k is the bulk viscosity 

coefficient, k is the coefficient of thermal conductivity, Cp is the spe- 

cific heat at constant pressure, u is the velocity, and T is the tempera- 
ture. In the effective one-dimensional flow just inside the holes, the large 
heat flux, -q , to the left Is accompanied by an effective viscous-compressive 
stress t. Unless Pr is zero or infinite, heat conduction in a compressible 
gas is in general always accompanied by a viscous stress . Recall also from 
kinetic theory (e.g., for a monatomic gas, see ref. 6, p. 1317 ) that the 

thermodynamic pressure is p = (l/3)p(cx 2 + 2c 2 2 ) and p - t = pc x 2 (where 
u + Cjl and c 2 are the components ofjnolecular velocity in the x direction 
and its normal, respectively, and ( ) is the average weighted by the dis- 
tribution function, so that t = (2/3)p(c 2 2 - c x 2 ) . The viscous stress r is 

different from zero if the averages of the different components of the random 

velocity are not equal. They evidently are not equal in the presence of a 
large heat-conduction flux, q c = pc 2 c x , which itself is different from zero 
because of the asymmetry of Cj_. (Here c 2 = C! 2 + 2c 2 2 . ) The asymmetry 
of Cx in general makes its average squared value different from c 2 2 . 
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Furthermore, the viscous stress just inside the holes cannot end abruptly, but 
must decay within a finite distance that is determined by the equations of 
motion. Therefore the viscous stress, and an accompanying heat flux, exist 
also at b. Thus, if there is a significant amount of heat absorbed at the 
wall, viscous effects will be present at b, and the region of their decay is 
the boundary shock wave. 

In the case of rapid vapor ablation, consider only a material that 
vaporizes from a liquid surface and has no chemical reactions after the vapor- 
ization. Assume the radiation is absorbed and the molten material is vapor- 
ized within a few molecular diameters at s . The high heat conduction to the 
left within that region is developing the conditions that will produce a 
viscous-compressive stress in the vapor state: t = (2/3)p(c 2 2 - c i 2 ) (tor a 

monatomic gas). Note that under mild conditions (small heat transfer, low- 
speed vaporization), the vaporization takes place reversibly, and the vapor is 
considered to be in equilibrium at the end of the phase change. Under more 
extreme conditions (high rate of radiation absorption and heat conduction just 
inside the surface, producing the high rate of vaporization) one expects 
mechanical (e.g., translational) nonequilibrium to be significant in the 
vaporization process. Translational nonequilibrium is described macroscopi- 
cally by the concept of viscosity . A large rate of heat transfer (to the left 
in fig. 1 or fig. 3) at s and within e is accompanied by a time lag in 
passing the translational energy from the axial mode (x direction) on to the 
degrees of freedom of the motion in the directions normal to the x direction. 
The lag makes c x 2 > c 2 2 ; that is, because of the lag, the molecules with 
motion mostly normal to the x direction are not moving (on the average) 
quite as rapidly as those in the x direction (relative to the mass velocity 
u) , and are therefore tending not to break their intermolecular bonds as rap- 
idly. The viscous stress t is a measure of this lag. (For polyatomic mol- 
ecules there will also be a lag associated with the rotational degrees of 
freedom, also affecting r.) 

The stress t develops rapidly during the vaporization because, in a 
pure liquid state, the molecules are effectively in contact and there is prob- 
ably very little lag in transferring the energy to modes normal to the axial 
motion; but the lag develops as the molecules acquire the additional energy 
that will free them from the constraints of' the intermolecular forces, and as 
they begin to separate from each other. By the time the phase transition is 
completed (at b ) , t he viscous stress t can be represented by 
r = (2/3)p(c 2 2 - c i 2 ) = p du/dx. This stress is accompanied by viscous dissi- 
pation of energy and an associated heat -conduction flux which, in the vapor 
state, can be represented by q c = pc 2 c x = -k dT/dx. 

In the above discussion, from the microscopic (molecular) point of view 
we see that there is a tendency, which is measured by t, for the vaporization 
process to be retarded. This has the effect of a compression opposing the 
expansion of the phase change, so t < 0. Refer to figure 3, on which the 


4 



sketch and curves are idealized and qualitative.* Let the mass velocity deriv- 
ative, du/dx, he separated into two terms (effects of two contributions) , 


du _ / du\ + / du\ 

dx phase W/visc 


where 


(du/dx) phase 


the contribution representing the tendency of the molecules at 
a station x to break their intermolecular bonds and sepa- 
rate from each other 


( du/dx) v ^ sc the retarding effect of viscosity on the velocity change 

(effect of mechanical, or translational, nonequilibrium); 
or a resistance to the separation of the molecules (associ- 
ated with t) 

At b, the term (du/dx) phase has vanished and the term ( du/dx) v i SC has the 
value appropriate to a gas, and is given by T/p,. The effect of ( du/dx) 
between s and b is a downward shift from the curve of (du/dx)phase to give 
the resulting mass velocity derivative, du/dx. In order for the vaporization 
process to be retarded by viscosity, t must be negative , so that the velocity 
derivative (proportional to acceleration), du/dx, during the phase change is 
decreased by the viscosity. This macroscopic reasoning is therefore com- 
pletely in line with the results of reasoning on the molecular level, above. 


The viscous effects that develop in the vaporization process cannot end 
abruptly, but must decay within a distance determined by the conservation 
equations. That is, although the intermolecular bonds are assumed to be over- 
come very quickly, the lag in communicating the translational energy to the 
lateral degrees of freedom may take much more time to accommodate. The region 
in which these mechanical nonequilibrium effects relax (i.e., viscous effects 
decay) with distance from their source is the boundary shock wave. 


It may be noted that, in order for to be negative as described 

above, the density must undergo a slight "undershoot” at the end of the vapor- 
ization process, and the mass velocity u undergoes a slight "overshoot. " 
Thus, near the end of the vaporization process, the curve of velocity versus 
distance (fig. 3 ) has a maximum point. This overshoot in mass velocity is 
actually exceedingly small, and the slope of the velocity curve at 
x = x-^ = 0 + € is extremely small in comparison to, say, Au/e (slope of 
dashed line on curve of u vs x) . The overshoot in velocity may be explained 
qualitatively by referring to figure 3 - From figure 3 one can see the over- 
shoot in u, or the maximum point on the u curve, as corresponding to the 


*0n figure 3 the curve of velocity, for example, is highly distorted. 
The change of u (or p) during the phase change, from s to b, is several 
orders of magnitude greater than the change from b to e. Also the distance 
from b to e is several orders of magnitude greater than the distance from 
s to b . 
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point where du/dx goes through zero. It is the point where (du/dx)p] aase 
has come hack down to the relatively very small value equal to -(du/dx) v j_ sc 
before vanishing. 


The arguments that have been presented in this section are not intended 
as proof that the postulated viscous effects will occur. Although one does 
not expect that the viscous effects will have vanished at b, absolute cer- 
tainty as to whether or not t-^ will be different from zero in a given case, 
so that a boundary shock wave will occur, would depend on either adequate 
experimental verification or on a rigorous analysis by nonequilibrium kinetic 
theory of liquids and gases for the phase-transition region. Only such an 
analysis could supply an exact value for . 


The values of the stress 


and the heat conduction q would be 


T b — — ^c^ 

determined by solving simultaneously all the equations pertinent to a specific 
problem, as is done in the example calculation in the last section of the 
analysis. However, the general theory of the boundary shock can be worked out 


in terms of 
section. 


and q , 
c b 


without specific application, as in the following 


GENERAL THEORY OF THE PLANE LAMINAR BOUNDARY SHOCK WAVE 
IN A PERFECT GAS WITH Pr = 1 


The general theory of the laminar boundary shock wave is formulated 
independently of the details of the mass-transfer process under the assump- 
tions: (l) the flow is close enough to mechanical equilibrium that the 

Navier-Stokes equations give a valid approximation to the description of the 
macroscopic thermodynamic and flow variables; and (2) the geometrical struc- 
ture of the surface is of small enough detail that the velocity vector can be 
considered to be essentially one-dimensional (normal to the wall) in the main 
portion of the region of interest. This is equivalent to the assumption 
e « & in the previous section. Then the conditions across the boundary 
shock wave and the structure can be calculated from the Navier-Stokes equa- 
tions and depend only on the values of the macroscopic parameters at the vapor 
boundary b . 

The exact equations for conservation of mass, momentum, and energy in 
one- dimensional steady flow in a nonaccelerating coordinate system are: 



= 0 

(5) 


II 

fclfl 

(6) 

eE ( e + 5 u2 ) 

= + uf) 

(7) 
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where x is the distance to the right of the boundary b, e is the internal 
energy per unit mass, and f is the sum of the "surface forces" in the x 
direction on an element of mass. In the gas, 

f = - p + T (8) 


It is often convenient to discuss the equations in terms of the specific 
enthalpy , 


h = e + p/p 


(9) 


rather than the internal energy. The equations of state for a thermally and 
calorically perfect gas are 


p = pRT, h = c p T , e = c v T 


( 10 ) 


where R, Cp, and c v are constants related by 


G p/ C v - 7 ? Cp - c v - R 


(id 


If heat transfer by conduction and radiation are considered, we can write for 
the heat flux in the x direction (positive to the right): 


1 + q r 


( 12 ) 


where q r is the radiative heat flux. Although radiation is included, it 
will be assumed that the gas is transparent to radiation and that radiation 
emitted by the gas in the boundary shock region can be neglected. Thus, all 
radiation is assumed to be emitted or absorbed outside the region from b to 
e, and q r is therefore constant across the boundary shock. In the Navier- 
Stokes formulation the viscous stress and the conduction heat flux are given 
by equations (l) . 

The boundary conditions to be used with equations (5)> (6), and (7) are: 
at x = x b = 0 + : u = T = T b , q c = q c ^ (13) 

as X -*• oo: — -<• 0, — -*• 0 (l4) 

dx dx 
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For simplicity, in the following sections, it will be assumed that Pr 
(eq. (3)) is equal to unity. This assumption is not overly restrictive since, 
as Liepmann and his coworkers point out, Pr ^ 1 for most real gases (ref. 6). 
The result, as has been shown for shock waves and as will be shown in a later 
section for the boundary shock (eq. (57) )> i s that 

q c = ru, (Pr = 1) ( 15 ) 


at every point in the boundary shock. With this result, it will be convenient 
to present further results in terms of 



( 1 / 2 )pb u b 3 ( 1 / 2 )p b u b 2 


( 16 ) 


Conditions Across a Boundary Shock for Pr = 1 

The conditions across a boundary shock wave are analogous to the 
conditions across an ordinary shock -wave, including the Rankine-Hugoniot rela- 
tion, the Prandtl relation, and other normal-shock relations (see, e.g., 
ref. 10, pp. 57-64) and are obtained in much the same manner. 


Integration of equations (5), 

(6) 

, and. (7) and use 

of the remainder of 

equations (l) through ( 15 ) give 





p b"b 

= P e u e 

( 17 a) 

P b + V* 

2 _ . 

T b = p e + p e U e 2 

(171) 

O u b 2 + ~ 
2 D 7- 

p b 

lp b 

1 2 7 k 

2 u e + 7-1 p e 

(17c) 


where subscript e denotes a value as^ x 00, that is, outside the boundary 
shock. Rote that the assumption of Pr = 1 has eliminated t and q from 
equation ( 17 c), so that the total enthalpy, h + (l/ 2 )u 2 , is constant. 

To derive a relation analogous to the Rankine-Hugoniot relation, combine 
( 17 a) and ( 17 b) to obtain 


P b Pe 
p e u e 


u e - u b 


P b U b 


(18) 
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Then multiply (l8) by u e + u^ to obtain an expression for u e 2 - u-]- 2 , which 
may be equated to another expression from ( 17 °) with the resulting relation 
analogous to the "Hugoniot relation" (cf. ref. 22, p. 6): 



2y /Pe £b \ 

7-1 V p e p b/ 


(19) 


This can he rearranged into a form analogous to what Liepmann and Roshko call 
the "Rankine-Hugoniot relation” (ref. 10, p. 64): 


i 7+1 _ _^b 

+ y- 1 p h ~ Pb 

P b 7+1 P e T 1 
l" 1 p b P b 


( 20 ) 


We may plot p e /p^ versus P^/Pgj which depends on (-r^/p^) = (l/^/M^C^ 

(see fig. 4). Later remarks are given concerning positive and negative values 
of this parameter. We may call this plot a "Boundary- Shock Hugoniot Diagram," 
by analogy to the Hugoniot Diagram for a simple shock wave. 

Next is developed a relation for boundary shock waves which is analogous 
to (and, in the case -> 0, reduces to) the Prandtl relation for shock 
waves . Denoting the speed of sound as 


a 



( 21 ) 


we can write equation (l8) as 


Up. 


- u c 



T b 

7 u b + Pb^ 


( 22 ) 


and ( 17 c) as: 



- (1±T) a* 2 = C 

2 \7-iy 


(23) 
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(Equation ( 23 ) may be regarded as the definition of a* and the constant C. 
In cases where u-^ > a^, a* is the actual value of a in the boundary shock 
wave where u = a.) Note from equation ( 23 ) that 


a b 2 = (7 - 1)C - -2“ u-b 2 ; a e 2 = (7 - 1 )c - -ai u e 2 


which may be substituted into ( 22 ) to obtain 


-b - u e = a* 2 ^ 


(24a) 


Equation (24a) may be further manipulated to give finally 


= a* 2 + -2Z 


u e u b “ a + 


7+1 V p e " p b 


(24b) 


For the case of the simple shock wave (t-^ = 0 and the ’’boundary” - 00 ) y the 
Prandtl relation, which determines that flow through a shock wave must go from 
either supersonic to subsonic or vice versa (see ref. 10 , p. 57 ) is recovered 
from equation (24b) . For the case of the boundary shock, such a restriction 
is not imposed. 

Other useful relationships across the boundary shock can now be developed 
as follows: Equations ( 23 ) and (24b) give two expressions for (a^/u ^) 2 which, 

when equated, give an equation that can be rearranged into the quadratic form 


where 



2B = “4 ( 1 + ~^~Z + k c h 

7+1 V 7Mb 2 


A = 


27 ( 2 zk + _L 


7+1 V 2 7 7Mb ; 




% 


(25) 


( 26 ) 
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The solution is 



( 27 ) 


which is thus determined for given 7, and (or given A and B) . To 

find the pressure ratio , combine equations (l7 a ) a^d (17b) into the form 



1 + r^ 2 





(28) 


The temperature ratio is obtained from the thermal equation of state in ( 10 ) 
and the energy relation (17c) as 


Te 

Tb 


1 + 


7“1 ,, 2 
2 “b 


1 - 



( 29 ) 


The Mach number downstream of the boundary shock, M e , can be found by writing 
equation (17b) in the form 


_1 

7 Mb 2 


+ 


1 


^ (hi. + A 

p b u fa 2 u b \/Me 2 ) 


from which 


2 

7 Me 




(30a) 


As an alternative, one could also use 


Me 


2 


u e 2 7 RT bMb 2 

7RT e u b 2 


into which equations 


(29) and (l7 a ) may be 


substituted to give 
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2 


2 


Me 


2 


\ P e/ ^ 


1 + ^ Mb 2 


1 - 



(30b) 


where p /p is given by (27). 
b e 

The conditions prescribed by the equations in this section are possible 
only if the associated entropy production is positive; thus the second lav of 
thermodynamics must supplement the conservation equations, (17)- It can be 
expressed as 


% 

A n -s = As - ApS = s e - s, — > 0 

1 ^ b ObVb - 


(31a) 


vhere s is the specific entropy, Aj_s is the entropy production per unit 
mass due to the irreversible processes (cf. refs. 28 and 29) within the bound- 
ary shock wave, and AeS is the change in entropy due to heat transfer (exter- 
nal process) across the boundary shock; or 


AqS 



7Mb 2c h c > 0 


(31b) 


The state variable, entropy, can be obtained for a perfect gas from 


s 


s ref “ C P ln iji 


T 


- R in 


ref 


• ref 


(32) 


With use of s^ as the reference value and s e as the value of interest, 
equation (32) leads to 


7"1 

R 


(s e 


- S h ) = Zn 


£e 

L P b 



(33) 


Finally, equations (31) and (33) 


give as the required condition 


7-1 

R 


A^s 



+ \ 7(7 - D\\ c > 0 


(34) 
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An interesting point here is that the dimensionless number on which the curves 
on the boundary- shock Hugoniot diagram (fig. 4) depend (in addition to 7 ) is 
just the ratio of the change of specific entropy of an element of mass due to 
transport of entropy (i.e., exclusive of entropy production within the mass 
element) across the boundary shock, to the gas constant R, that is, (cf . 
eqs. (31)) 



- \ ’• M b s °h c 


AgS 

R 


(35) 


For given 7 and Cjj the quantities [(7 - l)/R]Aj_s, [(7 - l)/R](s e - s-^), 
Pg/P^j ari< ^ ^e can now calculated and plotted versus 

using equations (25) through ( 3 * 0 . The results for > 0 are shown in 

figures 5 through 10 . The possibility of cases for which < 0 will be 

discussed later. Although the results shown on the figures were computed for 
7 = 5/3 * they are qualitatively the same for other values of 7 and are eas- 
ily computed from equations given herein. 

It is important to note first of all that the density ratio given by 
equation (27) and the other quantities given in figures 5 through 10 are 
double- valued because the quadratic equation (25) has two roots. The second 
law of thermodynamics (in the form of eq. ( 34 ) ) determines that only the minus 
sign in (27) is possible since, as seen in figure 4 , the roots corresponding 
to the plus sign, which are shown as dashed lines, are characterized by nega- 
tive entropy production. As a matter of interest, the total specific entropy 
change Is shown in figure 6 to point out that it is negative in some cases 
where the entropy production is positive as required by the second law of 
thermodynamics. Figure 6 is presented simply to illustrate this fact. It is 
interesting to compare figures 5 and 6. The difference between the total 
entropy change and the entropy production is, of course, the nondissipative 
part of the entropy change, that is, the transport of entropy due to heat con- 
duction (see. eq. (35)). 

From figures 7 through 10 we see that the case for which =0 is 

obviously singular. For the limiting case = 0 there is no boundary 

shock. If the efflux is subsonic (M^ < l) and -> 0 , the boundary shock 
simply vanishes. If the efflux is supersonic (M^ > l) and C^ c -> 0 , the 

results correspond to a simple shock wave with no boundary, that is, the bound- 
ary shock becomes a detached simple shock wave. (This latter case, the 
detached shock wave, is known to have been observed experimentally within the 
gas ejected supersonically from a porous blunt body (ref. 30 ). This has been 
referred to as "shocking down" of the injected gas.) 

Another important point to consider in calculating the possible conditions 
existing across a boundary shock wave is that the density ratio given by equa- 
tion (27) and, consequently, the pressure and temperature ratios and the Mach 
number M e are physically possible only if 
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(36a) 


B 2 > A 
or 





(36b) 


For 7 = 5/3; the curve of the minimum possible C^ c for P-^/pg to be real 

for specified is shown in figure 11. Note that this consideration also 

does not rule out negative values of (except at = l) , but does limit 

their magnitude. 

Results for positive have been discussed above , for which the heat 

is conducted back through the outflowing fluid, in the same direction as the 
heat conduction occurs in the physical situation existing in a simple shock 
wave. We may now also consider the possibility of negative values of 

for which heat conduction is in the opposite direction to that in a simple 
shock wave, that is, in the same direction as the gas flow. As an example cal- 
culations were made for = -0.1. The density and pressure ratios and the 

entropy production are shown in figure 12. As seen from figure 12, the Mach 
number at the boundary is limited to values less than 0.7 6 or greater than 
1.57 for this case. Both roots in equation (27) are shown in figure 12, and, 
as seen, the branches of the curves corresponding to the two roots are con- 
nected at = O.76 and 1.57 since they have common values there. The 
entropy production for the branch corresponding to the plus sign in equation 
(27) is always negative, and is negative for part of the branch corresponding 
to the minus sign. (Negative-entropy-production portions of the curves on 
fig. 12 are characterized by dashed lines.) But there do exist "low- speed 11 
and "high-speed" portions for which AqS > 0. For the low portion Aj_s/R 

for the branch corresponding to the minus sign in (27) is so small as to be 
imperceptible on figure 7 for = -0.1, but it is positive near - 0. 

No claim is made here regarding the physical existence of the cases for 
negative for the following reason. As pointed out by Hayes (ref. 22), 

existence of a gasdynamic discontinuity (which the boundary shock is for 
Re ^ = 00) depends on the stability of the phenomenon, both in the large and in 
the small, in addition to satisfaction of the conservation laws of mass, momen- 
tum, and energy, the equations of state, and the second law of thermodynamics. 
The cases with positive Cj lc correspond to simple shock-wave solutions which 

are known to be physically realistic, that is, to be stable and exist in real- 
ity. Although the cases with negative Cj lc do not correspond to solutions 

known to be physically realistic, and hence have not been proven to be stable, 
the possibility of their existence is not ruled out. In a given application, 
the actual value of would be determined from appropriate equations, as, 

for example, is done in a later section in the application to rapid vapor abla- 
tion. It should be pointed out here that, as Cft -+ 0 with > 1, the 

c 


Ik 



boundary shock detaches. Thus, at least for M 5 > 1, one may not expect the 
boundary shock wave to exist for < 0 in some applications, since required 

jumps in the flow properties may be accommodated by a detached shock. 

Equations are presented above for calculating conditions across a boundary 
shock wave for given values of 7 , M^, and C^. However, in a given problem, 

for example, one discussed in a later section, one may wish to calculate the 
conditions for given 7, , and P e /p^. The following equations are presented 

for convenience in that case. An expression for [1 + (1/2)0^] from equa- 
tion ( 28 ) may be substituted into ( 25 ) and ( 26 ) to obtain 



(37) 


where 


2B' 


2 £e 

(7-1)]^ 2 P h 


A’ 


1 + 


2 

( 7 -l)M b 2 




(38) 


which has the solution 


— = - B’ + -J (B* ) 2 + A* 
0 

M e 


(39a) 


(where only the positive square root is chosen in order that p /p be posi- 

2 b o 

tive, since B T is always positive). In many cases (B ! ) may be much greater 
than A, so that it may be difficult to obtain an accurate number from (39 a ) • 
In those cases the root of (37) is computed more easily from 



A T 

B' + n/ (B* ) 2 + A' 


(39b) 


Calculation of Cv > Te/Th> an( ^- Me then follows from equations (28), ( 29 ), 

and ( 30 ). c 

In later work it will also be found convenient to make calculations for 
given values of 7 , and Mg. For that purpose the following equations are 

given. Equation (30b) can be solved for (P b /P e ) as 
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Mg 2 

V 


7-1 *, 2 
2 


1 + ^ % 


7-1 2 

1 + — M a 2 


Equation (29) becomes 



7-1 Mb 2 - Mg 2 
2^ l+r-lMe 2 


(40) 


(4l) 


Since the density ratio is known from (4-0), the pressure ratio can he computed 
from 


Pg = P e u e 2 7Mfa 2 _ p b Mb 2 

P h 7M e 2 P^ 2 P e Me 2 





(42a) 


(42b) 


and (1/2)0^ can be found from equation (28), with use of (4o) and (42a), as 



Structure of a Plane Boundary Shock 

In order to calculate the variation of the flow variables through the 
boundary shock wave, one may begin with the integrated forms of equations (5)j> 
(6), and (7) and with the boundary conditions (13) and (l4) . (Conditions (13) 
are actually just definitions of the symbols for the boundary values. To 
these could be added: p = t = t^, p = p^, etc., at x = x^ = 0 + .) With 

use of the appropriate auxiliary expressions, including the equations of state, 
(10), the differential equations and boundary conditions may be put in terms 
of only two variables, u and T. Then it is convenient to make the equations 
dimensionless by defining and substituting the dependent variables: 
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ii 


1 ia 11 


HUB! ■ II 


u = u/u b , T = T/Tb 




and a new independent variable 




5 


dx 

£ 


( 45 ) 


r» is con- 
P 


With the assumption Pr = const., £ is proportional to k, since 

stant, but the variation of p need not be further specified if the results 
are left in terms of the independent variable We now have the equations: 


du _ _ _ 
d| 


7% 2 u 


p hV 


- 1 




(46) 


d 

f 1 T ,1 ^ 2 ~| 


r t , 1 u 2 ] - 

r 1 + i] 



d| 

_Pr (7-DMb 2 2 J 


_( 7 -l)Mb 2 2 J 

_(y-D % 2 - 

Vb 3 

p b u b7 


(47) 


and the boundary conditions 


t 




0 


(48) 


In the special case considered here for which Pr = 1, it is convenient to 
define 


H = 


CpT + | u 2 
u b 2 


(7-1)1% 


1 _2 
+ — U 
2 2 


(49) 


Note, then, that 


dH 

cL| 


, dT ^ du 
k _ ~ q c + UT 

p b U b 3 + ~ p bV 


(50) 


17 



Further define: 


so that (47) is simply 


6 = H - Ik 


ae _ a _ ( de 


d| 


At, 


where, by definition. 


at t = 0, 0 = 0 


and where, from the boundary condition (48), 


d0 

as | -» 00 — 0 

B J d| 


The only possible solution is, of course. 


0=0,^ 

at 


AS 

as 


Hence, equations (49), (50), and (51) give 


T = 1 + ^ ^ [l - ti 2 


and 


q = ut 
c 


from which also equation (l6) stems. Equation (46) then becomes. 


— = au - b + - 
d| u 


where, by definition. 


(51) 


(52) 


(53) 


(54) 


(55) 


(56) 


(57) 


(58) 
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at £ =0; 11 = 1 


( 59 ) 


and the boundary condition is: 


I 



0 


(6o) 


and where 


a 


7+1 

2 7 


b 



7-1 + _1 


(61) 


(There should be no confusion between this constant, a, and the speed of sound, 
a, used earlier.) The solution for d 2 > 0 (see eq. (36a)), where 

d = s/b 2 - 4ac = 2a n/b 2 - A (62) 


is 


,- 2a i _ 


a - b + c 


- bu + c, 


au 2 


/ 2au - b + d . 2a - b - dV 3 ^ 
\2au - b - d 2a - b + d / 


( 63 ) 


from which , given 7, M^, and Ch c , one can calculate u versus £ (or, given 
just the two parameters b/a = 2B and c/a = A, one can find u versus a£; 
compare_(25) and ( 58 ) for the case u = ue ; du/d£ = 0) . Then equation ( 56 ) 
gives T versus £ , and finally one can find 


and 



( 64 ) 


versus Results for a typical case (7 = 5 / 3 , = 1, = 1/2) are shown 

plotted in figure 13 . 

The variation of £ with x, and hence of the flow variables with x, is 
easily found if ft is given as a function of T. This variation has been 
discussed by many authors in connection with the shock- wave problem. The var- 
iation of p as a gas property, for which reference can be made to the 
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literature , is not considered here, as the intent was simply to show how the 
houndary-shock-wave structure could he calculated, in addition to obtaining 
the result (57) needed in above developments. 

APPROXIMATE SPECIAL APPLICATION OF THE 
GENERAL THEORY TO RAPID VAPOR ABLATION 


The purpose of this section is to make a calculation of the boundary 
shock wave near the forward face of a blunt body having the properties of a 
stone meteor entering the earth 1 s atmosphere under conditions for which the 
boundary-shock phenomenon might be expected to occur. 

In an effort to eliminate all unnecessary complications, a number of 
idealizations and approximations are used in the present application. Con- 
sider the same assumptions used in the development of the theory above. The 
applicability of the properties of a perfect gas with Pr = 1 in the vapor 
from a stone meteor will be discussed later. The assumption of steady flow 
is here one of quasi-steady flow, for which time derivatives are neglected in 
comparison to space derivatives, with forces in the body, or wall mass, 
replaced by accelerations of the body. The ablating surface is considered to 
be receding back into the liquid material with constant relative speed u s , so 
that in the coordinate system of figure l(c), the surface x = 0 is station- 
ary and the wall material is moving toward the surface with constant velocity 
u s . Assume that very little heat radiation is emitted or absorbed within the 
short distance to the right of b in figure l(c), but that all radiation is 
absorbed at the surface, within the very small distance € in figure l(c) . 
Assume also that no chemical or phase changes take place to the right of x-^ 
and that only vapor is blown off. 

With very large Re^, the region of possible rapid variations under 
consideration in the flow near the surface is very thin, so that consideration 
of one-dimensional flow out of a plane surface gives a valid approximation and 
conditions (l7) are the appropriate first-order outer boundary conditions. 

The one-dimensional equations of motion which apply to the material regardless 
of its phase, and hence which can be used to relate conditions to the left of 
x = 0 (in the unvaporized state) to those in the flowing- vapor state are given 
by (5) , (6), and (7)> or in their integrated form: 


pu = m = constant 
pu 2 - f = constant 
1 ^ 

pue + — pu 05 + q - uf = constant 


( 65 ) 


Equations (l) through (15) apply for x > 0. The region from s to b in fig- 
ure 1(c) is treated as being infinitely thin, so there are jump discontinui- 
ties in velocity and density at x = 0 due to the phase change. To the left 
of x = 0, the compressive stress is -f, the equation of state is 


20 



( 66 ) 


p = p = constant, (x < 0) 
s 


and the heat flux is 


<1 = «L C 



(x < 0) 


(67) 


where k is appropriate to the state of the material at the particular loca- 
tion x. With subscripts a, f, and s denoting values at the respective 
locations: arbitrary point inside the wall at the cold solid state (x a ) , 

fusion interface (xp ) , and surface just inside radiation-absorption and vapor- 
ization interface (x s = 0~), one can write for the internal energy in x < 0: 


e 

e 


e a = c soZ (T - T a ), (x < x f ) 

e a = c soz( T f " T a) + L f + c Uq (T - T f ) , (x f < x < x s = 0") 


(68) 


where c sol a nd c Ziq are s P ec i f i c heats in the solid and liquid, taken to 

be constant, and is the latent heat of fusion. 

For x > 0 the solution to equations ( 65 ) is given in the above section 
on the general theory of the boundary shock wave . 


For x < 0 equations ( 65 ) 9 ( 66 ), and ( 67 ) give 

-\ 





which simply determines the temperature distribution inside the wall and, 
hence, relates the heat conduction and temperature at the surface (x = 0~) to 
the heat conduction and temperature at point a inside the wall. By defin- 
ition, x a is in the cold interior of the wall, where T = T a = constant and 
dT/dx = 0, so the term k a (dT/dx) a in ( 69 ) is zero. 

To relate the conditions across the surface at x = 0, equations ( 65 ) 
take the form: 
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P b u b = P s u s = m 


P^ u b 2 + Pk " T b = P s Us2 “ f s 


( 70 a) 

( 70 b) 


P b u b e b + \ P b u b 3 + q Cfe + q r + u b (p b - r b ) = p s u s e s + ± p g U s 3 - (k ^ 


^ — u s f s 

( 70 c) 


Noting from equations (70a) and (70b) that 


fs = -P b V U b - U S ) - P b + 


(71) 


and substituting this into (70c), we find 


' q r = (7 f") + % U b(% ~ e s) + \ P^K " u s) 2 + (P b ~ T b )^ b - u s ) 


(72) 


The internal energy change across x = 0 is approximately 


~b e s ~ 


(73) 


for vaporization either at low rates or at high rates, where is the 

latent heat of vaporization for low rates (see appendix B) . It is convenient 
to define (using eqs. (68) and (73)): 


5 ab 


= e b “ e a + RT ‘ 


(7^a) 


- c soz( T f ~ T a) + L f + c ziq( T b ~ T f) + W 


(7^b) 


which is the heat that must be transferred to a unit mass to raise it from the 
cold state at temperature T a to the vapor state at temperature T^, includ- 
ing all phase changes. (Note that T^ has been taken to be equal to T s ; 
see appendix B.) Then the terms 



+ P b V e b 


(75) 


I 
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in equation ( 72 ) can be replaced, using ( 69 ), ^y 


V b ( S K - ET b) (76) 

Equation ( 72 ) is further simplified by using 

u b ■ U S = u b( X - P b / p S ) » u b ( 77 ^ 

to obtain 

“4 r = p b U b (j>ab + | u b 2 ) + % “ T b u b ( 7Q ) 


For Pr = 1, it was shown above that the solution for the vapor flow is possi- 
ble only if q c = tu (eq. ( 57 ) )> with which we have finally 


= p b U b 



,(&-!) 


(79) 


where £ , given by equation (74b) , is a linear function of T b . 

a d 

For given values of the incident radiative heat flux (-q ) and the 
external pressure (p e ), the above equations (including the equations for the 
boundary shock in the previous section) are sufficient to determine the 
boundary- shock solution in terms of two unknown quantities, say p^ and p^u^. 

Some approximate equations describing the vaporization process may be used to 
relate p^ and p u-^ to the other parameters of the problem, and thus to 

complete the determination, as described in the remaining paragraphs of this 
section . 

The rate of vaporization, or the mass flux, depends on the velocities of 
the molecules leaving the surface during the phase change, and hence on the 
temperature at vaporization. For present purposes, let us accept the roughly 
approximate relation used by Opik (ref. 31* P* 24): 

(pu) b = | P b “ x (80a) 


where p^ is the vapor density at saturation pressure and (l/ 2 )v x is the 
average component of molecular velocity in the +x direction. For typical 
meteoritic stone, this reduces (ref. 31 * P • l 6 l) to: 
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K 

i 



1 gm cm ^sec" 


1 dyne/ cm 2 / \% / 


kY /2 


With use of the equation of state (10), equation (80b) gives directly 


i/Wb v 


= 3-08 x 10' 


R/y 

1 cm 2 sec“ 2 (°K) -1 


(80c) 


For meteoritic stone, Opik gives (ref. 3U P* l6o) values for the mean atomic 
weight as 23 and the mean molecular weight of vapors as 50. Hence , 

R = Ru/( 50 gm/mole) = 1.66 X 10 6 cm 2 /sec 2 °K, where is the universal gas 

constant. Since the mean molecular weight is about twice the mean atomic 
weight, we assume the gas mixture can be approximated as a perfect diatomic 
gas with 7 = 7/5 and c p = 7^/(7 “ l) = (7/2)R = 5-8 X 10 6 cm 2 /sec 2 °K. 
(These values seem to be realistic, since Cp for the vapor is then about 
half the value of the liquid specific heat given by Opik as 1.1 x 10 T . By way 
of comparison, the specific heat, Cp, of water vapor is about one-half the 
specific heat of liquid water.) The result from equation (80c) is then 
^ = 0.333. 

The vapor pressure, p^, and the boiling temperature, Tb, may be related 
approximately by Clapeyron*s equation: 

Lv/ 1 l\ / n 1 \ 


R VT r , 


ref x b 


Clapeyron's equation is derived by assuming that the phase transition takes 
place reversibly at constant temperature and pressure while the entropy and 
specific volume (or density) change (ref. 32, p. 317; see also appendix B) . 
Strictly speaking, therefore, Clapeyron’s equation would apply only to vapori- 
zation at low rates. The use of Clapeyron f s equation as a reasonable gross 
approximation relating p, to Tp in the present application with higher 
vaporization rates is justified because of: (l) the exponential (rapidly 

varying) dependence of Pb on T^ (see appendix B) , and (2) the fact that if 
one makes a gross intuitive modification of Clapeyron's equation to account 
for the viscous effects and the variation of pressure during the vaporization 
process, the results of the calculations are not significantly different from 
those obtained using Clapeyron's equation (see appendix B and numerical 
results given below) . 

For the application to the high-speed meteor composed of the typical 
meteoritic-stone material, Opik (ref. 31, P- l6l) gives the mean boiling point 
as 2960° K (at 760 mm Hg or 1.013 X 10 6 dynes/cm 2 ) which can be used as the 
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reference value in equation (8l) . The value for R is estimated above 9 and 
L v can be taken from Opik T s data (given below). For reasons that 'will become 
apparent^ it is expedient to re-write equation (79) as a relationship between 
q r and T^ in terms of known quantities; using equation (8l) in the form 


(82) 


(83) 


The result is 


and the equations 




1 dyne / cm* 


= exp 


(26.13 - .35,500 °K 


P b = P b /RT b and u b = ^ 


-1, 


1 erg cm 2 sec 1 


expf 26.13 - 



36,500 


f W T b) 


1 cm 2 sec~ a 


RT> 


+ | 7Mf 


(84) 


where £ = (T b ) is given by equation (74b) and where, from Opik (ref. 31, 

pp. 6l and l6l), values of the required constants in (74b) for this case are: 


L f 

c so2 

T a 


2.65 X 10 9 erg/gm, 
8.95 x 10 s erg/gm °K, 


200° K 


W 

c Ziq. 

T-p 


6.05 X 10 10 erg/gm 
1.1 x 10 7 erg/gm °K 
18 00° K 


If the radiative heat flux to the wall causing the ablation, -q r , and 
the external pressure, p , are specified, then the solution of the equations 
for the boundary shock wave is determined and can be calculated by iteration 
on equation (84). Thus, a trial value of Tfc, may be chosen and iterated upon 
until the proper value for q r is obtained. Then p b is determined by equa- 
tion (82); the ratio P e /P b can be computed; and equations (39), ( 28 ); (29); 
and (30) can then be used successively to find p e /p^; Cp c ; Te/'Pb; an( ^- Me* 

The radiative heat flux, -q ; is a function of conditions of a particular 
problem. It depends ; for example; on a body radius; speed; and altitude. 
Representative values of -q r and p e used are; respectively; 17; 300 watts/cm 2 
and 7.55 X 1G 5 dynes/cm 2 ; found from a simplified approximate calculation of 
flow in the forward region of a spherical body of 4.62 m radius moving at 
15 km/sec at an altitude of 60 km in the earth* s atmosphere. The calculations 
as outlined herein; using the numbers given above corresponding to stone 
meteors; lead to: = 0.333; T^ = 2720° 


K; p^ = 3.345 X 10 5 dynes/ 


cm 2 ; 
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P b = 7.^1 X 10 5 gm/ cm 3 , ut, = 2.65 X 10 4 cm/sec, Pg/l^ = 2.257, P e /p b = 2 « 21 5, 
Mg = 0.1^9, and = 13.772. If the "modified Clapeyron's equation" is used 

(appendix B) , the above numbers are the same except for replacement of corre- 
sponding values by: -q r = 17^120 watts/cm 2 , r-^ = ^-.5 8 T-^ = 2,880° K, 

= 7-00 X 10” 5 gm/ cm 3 ; and u b = 2J3 x 10 4 cm/sec. 


CONCLUDING REMARKS 


An attempt has been made to present an introduction to a possible 
phenomenon defined herein as a boundary shock wave, a thin layer of gas 
through which viscous effects decay in very high-speed mass transfer accompa- 
nied by large heat transfer. The general theory of laminar flow through a 
plane boundary shock was developed. Results from the general theory include 
calculations of the conditions existing across a boundary shock, and of the 
internal boundary-shock structure, for given conditions at the boundary. The 
solution of the governing equations was reduced to closed-form expressions for 
the quantities of interest under the restrictions of a perfect gas flowing at 
constant total enthalpy. Numerical results were obtained and plotted for the 
case where the specific heat ratio, 7, is 5/3 • 

This theory -was then applied to rapid vapor ablation from a body with 
flight conditions appropriate to a possible meteoric fireball in the earth's 
atmosphere . 


Ames Research Center 

National Aeronautics and Space Administration 
Moffett Field, Calif., Sept. 22, 1965 
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APPENDIX A 


A, B 
a 


a* 


C 


C h 


c 


C l; C 2 

C V; c p 
e 

f 

H 

h 

k 

L 

Lf ^ Lv 

lv 

M 

m 

Pr 

P 

a 

^c 

^r 


PRINCIPAL NOTATION 


constants defined by equations (26) 
speed of sound (eq. (21)) 

critical speed of sound defined by equation (23) 
constant defined by equation (23) 


conduction heat-transfer coefficient. 


1 

2 



3 


components of random molecular velocity in the direction of x, and 
perpendicular to x, respectively 

specific heats at constant volume and pressure, respectively 

specific internal energy 

stress, or "surface force" per unit area 
dimensionless total specific enthalpy (eq, ( 49 )) 
specific enthalpy (eq, (9)) 
thermal conductivity 

characteristic length, especially width of vapor region 

latent heats of fusion and vaporization 

viscous length, q^/p^u^ 

Mach number, — 

7 a 

mass flux (eq. (65)) 

_ £ c p 

Prandtl number based on ll, — - — 

k 

thermodynamic gas pressure 

heat flux (to the right in fig. 1 (c)) 

conduction heat flux (eq. (l)) 
radiative heat flux 
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gas constant (eq. (10)) 

Reynolds numbers defined in figure 1 

body radius 
specific entropy 
temperature 

dimensionless temperature (eqs. (44)) 
time 

velocity in the x direction 
dimensionless velocity (eqs. (44)) 

space coordinate, positive to the right in figure 1(c) 

change in specific entropy due to heat transfer (external 
process) across a boundary shock (eq. (35)) 

change in specific entropy due to irreversible (internal) 
processes, that is, entropy production, across a boundary 
shock (eqs. (3l)) 

total change in specific entropy across a boundary shock, 
s e — s b (eq. (33) for a perfect gas) 

C P 

ratio of specific heats, — 

c v 

heat that must be transferred to a unit mass of material to 
raise it from the cold state at temperature T a to the 
gaseous state at temperature (eqs. (74)) 

defined by equation ( 51 ) 

defined by equation (4) 

dimensionless independent variable defined by equation (45) 
mass density per unit volume 
defined by equation (l) 



Subscripts 


a 

b 

e 


s 

1 , 2 


* 


value inside the wall (x < 0) where k — 0 

dx 

value at the boundary in the gas (x = 0 + ) 

value in the gas outside the boundary shock wave (x co) (or outside 
the boundary layer in fig. l) 

value at the surface , inside the wall (x = 0“) 

values upstream and downstream, respectively, of a simple shock wave 

Superscript 

condition where u = a 
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APPENDIX B 


APPROXIMATE THERMODYNAMIC EQUATIONS FOR RAPID VAPORIZATION 


Two different approximate formulations for the vaporization process are 
considered. The first is well known, but, for quantitative accuracy, its use 
might be limited to fairly low rates of vaporization. However, the results 
serve as a foundation for the proposed formulas applying to more rapid 
vaporization. 

The first law of thermodynamics for an element of mass undergoing a phase 
change is 


T ds = de + p dv 


(Bl) 

& 

> 

i 

£ 

ii 


(B2) 

= c v dT + T 

!E) dv 

ay v 

(B3) 


where v = l/p is the volume per unit mass. Denote the initial state 
(liquid) by subscript i and the final state (vapor) by subscript f . 

1. First approximate method. - 

If the phase change takes place at a low rate, the process is nearly 
reversible, isothermal, and isobaric. The heat transferred reversibly to the 
element, per unit mass, or the latent heat of vaporization, is, from 
equation (B2), 


Lv - T(s^ - s i ) - h^ - h.j_ 
= e f “ e i + P( v f _ v i) 
« e f " e i + P v f 
= + RT 

Since p and T are constant during the process, 

P = P(T) 

so that 



(B4a) 


(B4b) 


(B5) 


(B 6) 
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Integration of equation (B3) for the process at constant 


T gives then 


L v = T — (v f - v ± ) Rs T 


dp 


dT 


dT 


(B7) 


Use of the perfect gas equation of state: 


Pf v f = RT f (B8) 

in (B7) then gives 

^ = h. M (b 9 ) 

P R T 2 


Equation (B9) integrates (for constant L v ) to give the well-known Clapeyron's 
equation for the pressure to he overcome by the molecules escaping the liquid 
(the saturation vapor pressure) as a function of the vaporization (boiling) 
temperature: 


P 


f 


P 

ref 


Ly ( 1 „ 1_\ 

e R vref T f / 


(BIO) 


For a vaporization process rapid enough that it is not reversible and the 
pressure is not constant, Clapeyron’s equation should still provide a reason- 
able first approximation for p^ as a function of for the following rea- 
son. The variation of p f with is exponential, so the variation of p f 

due to changes of p during the process and due to irreversibility will be 
negligible in comparison to the dependence on T if those effects are not 
large . 


2. Second approximate method.- 

For higher rates of vaporization the process is not reversible because of 
viscous dissipation, and not isobaric because of significant viscous stress, 
which makes the normal stress on a fluid element different from purely thermo- 
dynamic pressure. The following estimate of the relationships among the ther- 
modynamic variables for the vaporization process is proposed as an improved 
approximation in that case. 

The steady- flow momentum equation is 

p + pu 2 - t - constant (Bll) 
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For the first approximation above it was implicitly assumed that both pu 2 
and ~t were negligible, so that p was constant during the process of the 
phase change. For this improved approximation, assume only that 


pu s « p - T 


(B12) 


(with this condition remaining as a limitation on the use of the final results). 
Then, during the vaporization, we have, approximately, 


p - t = p f - T f (B13) 

For rapid vaporization, the total heat supplied to a unit mass during 
vaporization could be denoted, using equation (B2), as 


Iv* = T ( s f ~ s i) “ kf " 



(Bl4) 


Although Lv* may be quite different from the value for slow vaporiza- 

tion rates, equation (B4b) may still be retained as a valid approximate 
energy relation for high rates. The reasoning is as follows: In the process 

of vaporization from a liquid surface, the individual molecules escaping the 
intermolecular bonds (close-range forces) acquire a certain amount of potential 
energy , which, in the macroscopic sense, is thought of as increased internal 
energy of the fluid. This change in internal energy of the fluid, or the 
value of the potential energy added to the individual molecules, is assumed 
not to be dependent on how fast the interface between the vapor and the liquid 
recedes into the liquid, that is, on the rate of vaporization. Thus the same 
internal-energy change is required for rapid vaporization of the material as 
under conditions of low-speed vaporization at the same temperature. Since 
« p^, the enthalpy change is approximately 

h f - h ± = e f - e i + RT (B15) 

Thus, for the same temperature and same (e^ - e ±) > the enthalpy change is the 
same for high rates as for low rates; that is, 

h^ — hq = e^ — eq + RT — (Bl6) 

An additional consequence of the above arguments is that the temperature 
remains nearly constant during vaporization at a high rate even though the 
pressure is changing and viscous effects are important. Since hq = h^(T) 

and hf = hf(T) and since hq - hq is the same for high rates as for low 
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rates at a given initial enthalpy, , and Tf = for low rates, it follows 

that T f = Tj_ for high rates of vaporization as well. 

In order to understand more completely the relationship between L v and 
L v * for rapid vaporization, consider the following. Equations (Bl4) and 
(Bl6) , along with (B13) , show that 


fPf 

Lv - By* - v dp 
J P ± 



(BIT) 


The quantity L v * is different from the latent heat for low rates, L v , 
because of the variation of p during the process at high rates , or ; equiva- 
lently, because of the significant variation of t. The relationship between 

P T f 

/ v dT in equation (B 17 ) and the irreversible heating due to viscous dis- 
^ o 

sipation in the process can be found by first writing, for the heat addition, 
per unit mass, to the mass element considered in equations (Bl) to (B3): 


T ds = T d e s + T dqs 


(bis) 


where 


T d P s = 


-dq -dq c - dq r 


pu 


pu 


(B19) 


is the heat transferred to the differential element per unit mass during the 
vaporization process. 


q = q r + q c (B20) 

is the heat flux (positive in the direction of flow) composed of the flux of 
radiation being absorbed during the process and flux of heat conduction, and 
where 


T dqs 


t du 
pu 


(B21) 


is the heat added by viscous dissipation per unit mass of the differential 
element. Now, from equation (B17)> 


Ly - Ly 


* 




u f 

_ r % du 

pu J., . pu 
U 1 


(B22) 
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where pu is constant for the vaporization process , by conservation of mass. 
With equation (B21), then, 

t p( f ) 

Ly - Ijv* = ~ _ / T d iS (B23) 

P f J (i) 

Thus the difference between Ly and Ly* is due to the work term 

TfUf/pu = T p/pp as well as to the irreversible heating by viscous dissipation. 

Even if the irreversible part of T(sp - Sp) in (Bl4) is much less than the 
reversible part (which is true because of large q r being absorbed in the 
vaporization process), Ly* can be significantly different from Ly because 
the reversible part of T(sf - s±) is different due to the pressure change. 
However, as pointed out above, hp - hp = L v remains the same at a high vapor- 
ization rate as at a low rate. 

The vapor pressure, p f , in Clapeyron T s equation, (BIO) , is interpreted as 
the "pressure to be overcome 11 by the molecules escaping the liquid. Since for 
rapid vaporization there is an additional stress, -Tp, to be overcome, let us 

assume that Pp, in Clapeyron T s equation can be replaced by p f - Tp to 

obtain a rough approximate relationship between the final total compressive 
stress of the process, p f - Tp, and the temperature, Tp = Tp: 

„ _ r i A 

f Tf _ e R VT r ef Tf J (B24) 

P re f 

Recall from kinetic theory (for a monatomic material) that if the material is 
in equilibrium in the final state, as assumed in deriving Clapeyron T s equation, 
then, 






(see the section "Causes of Boundary Shock Waves" for terminology). Also from 
kinetic theory in the case of mechanical nonequilibrium (i.e., significant 
viscosity) , 


(P - T) f = (pC! 2 ) f 

Hence, if (pci 2 )p is regarded as the quantity limiting the outflow of new 
vapor molecules (equal to the pressure in the equilibrium case), then it is 
natural to write as an approximation the above modified Clapeyron equation, 
(B24), for the nonequilibrium (viscous) case. Since the relationship between 
Pp and Tp now involves Tp, for a given application Tp would have to be 
determined from the coupling with the equations for the flow outside the imme- 
diate region of the phase change. 



I 


Although the intuitive result, equation (B24), was not derived rigorously, 
it is consistent in that, as Tf 0, Clapeyron's equation is recovered. The 
use of this "modified Clapeyron f s equation,” (B24), is justified by the fact 
that there must be some such relationship characterizing the properties of the 
material during the rapid vaporization, and the results obtained with its use 
must be assumed to be qualitatively, rather than quantitatively, correct. 

Also, equation (B24) is useful in determining when Clapeyron's equation, (BIO), 
is valid quantitatively. If under given conditions both (BIO) and (B24) give 
nearly the same quantitative results, it can be assumed that those results are 
accurate. 
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Figure 1.- Analogy of boundary shock wave with boundary layer and shock wave, 
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Figure 3*- Development of viscous stress during rapid vaporization. 
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Figure 13.- Boundary-shock structure: variation of pressure, density 

temperature, and velocity within a boundary shock for 7 = 5/ 

Mb = 1, C hc = 1/2, Pr = 1. 
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